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Abstract 



The inner crust of a neutron star is a superfluid and inhomogeneous sys- 
tem, consisting of a lattice of nuclei immersed in a sea of neutrons. We 
perform a quantum calculation of the associated pairing gap and compare it 
to the results one obtains in the Local Density Approximation (LDA). It is 
found that the LDA overestimates the spatial dependence of the gap, and 
leads to a specific heat of the system which is too large at low temperatures, 
as compared with the quantal result. This is caused by the neglect of proxim- 
ity effects and the delocalized character of the single-particle wavefunctions 
close to the Fermi energy. It is possible to introduce an alternative, simple 
semiclassical approximation of the pairing gap which leads to a specific heat 
that is in good agreement with the quantum calculation. 

PACS numbers: 21.60.-n, 26.60. +c, 97.60. Jd 
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I. INTRODUCTION 



Above a density of about about 3 x 10^^ g cm~^ (corresponding to 1.8 xlO~'^ n fm~^) 
matter in neutron stars consists of a matrix of nuclei immersed in a sea of neutrons and 
an approximately uniform sea of electrons. Such a configuration persists up to roughly half 
nuclear saturation density, and constitutes the so-called "inner crust" of a neutron star. 
At low temperatures, this system is superfluid with a positive Fermi energy ep-, and one is 
confronted with the task of estimating the pairing gap A of nuclei immersed in the neutron 
liquid [Q . With the exception of ref. 0, the pairing gap in this inhomogeneous medium 
has been calculated assuming that locally it coincides with the value of the pairing gap in 
bulk matter at the corresponding density (cf. e.g. ref. and refs. therein). 

A complete calculation of pairing should include the induced interaction, taking into 
account polarization effects in the medium (cf. for example 0). In the present paper we 
will restrict ourselves to the mean field level, and we shall compare the results of a local 
density approximation (LDA) calculation of the pairing gap to the results obtained from 
a Hartree-Fock-Bogoliubov (HFB) calculation. We find that the LDA leads to a spatial 
variation of the gap near the surface of a nucleus which is stronger than that obtained in the 
HFB calculation. This is because a Cooper pair in the quantum (HFB) calculation samples 
the pairing interaction over a relatively large distance, which is of the order of the coherence 
length (cf. Section HI). As a result, the pairing interaction in the neutron liquid influences 
the gap in nuclei, and vice versa, so that the spatial dependence of the gap is smeared out 
near the nuclear surface. 

The pairing gap plays an important role in the thermal evolution of a neutron star, 
because it influences in an essential way the specific heat in the inner crust. We will show 
that the LDA overestimates the effect of the presence of nuclei on the specific heat by a large 
amount (cf. Section IV). In the LDA, the specific heat is obtained from an integral in phase 
space, and it is possible to single out two contributions to the specific heat, one from the 
interior of the nucleus and one from the outer neutron liquid. Inside the nucleus the local 
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pairing gap is small, and this enhances the value of the specific heat. On the other hand, 
in the quantal calculation the specific heat is determined from the value of the pairing gap 
of the states close to the Fermi surface, which is obtained from an integral over the whole 
Wigner cell, and is much less sensitive to the presence of the nucleus. 
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II. PAIRING FIELD IN COORDINATE SPACE 



We consider a spherically symmetric system which is governed by a single-particle Hamil- 
tonian h and a two-body interaction v. The Hamiltonian h = T + U consists of the kinetic 
energy and a mean-field U which we do not calculate self-consistently but parametrize as 
described later on. The explicit effect of the two-body interaction on the ground state of 
the system is therefore included in the HFB equations only via the pairing field. 

The system is placed into a spherical box with a nucleus at the origin. The radius of the 
box is equal to that of the Wigner-Seitz cell of the lattice (cf. Fig. 1). Nuclei belonging to 
different cells are so far away from each other, that one can solve for a single cell separately. 
The single-particle wavef unctions, which are eigenfunctions to h are therefore required to 
vanish at the boundary of the cell. They are denoted by (pnijmir) = r\nljm >, while enij are 
the associated eigenvalues. 

The HFB equations 0] are solved in matrix form using this basis. We shall, in particular, 
make use of the two-particle wavefunctions 

< rir2\nn'lj, J = > = 0n«i(ri)0n'«i(r2) Yl r, ■ , . < hlljm >< faj/j - m >, (1) 

where J is the total angular momentum of the system and (pnij is a radial wavefunction. 
The normalization of these two-particle wavefunctions is 2/(2j -|- 1). 

In the case of spherical symmetry, the HFB equations are diagonal in the quantum 
numbers (Ij), having the form 

(e„ij - (^F)U^lj + Yl ^nn'ljV^'lj = EqijU^^j, (2) 
n' 

- {(-nlj - ^F)V^lj + ^nn'ljU^'lj = EgijV^ij, (3) 
n' 

where ep is the Fermi energy. The eigenvalues Eqij of these equations are the energies of the 
quasiparticle states, while the U and V eigenvectors are the occupation amplitudes on the 
single-particle states. 
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The matrix elements of the pairing field, An^/jj, can be expressed as 

^nn'ij = - < nn'lj; J = 0\v\^ >, (4) 

where 

1$ >= E ^^^(E U'ni.Kl,^) Inn'lj; J = 0> (5) 

nn'lj Q 

is the pairing density, which is related to the abnormal density as discussed in ref. 0. 

The HFB equations are solved self-consistently. Each set of solutions {U, V) determines 
a new pair density $ and new matrix elements of the pairing field, which are then used 
in the next iteration. This scheme is continued until a converged set of solutions has been 
obtained. The two particles in Eq. can be in states with a different number of nodes, 
while in the BCS approximation they are assumed to be in time reversed states (n = n'). 
The condition n = n' is in general too restrictive for an accurate description of localization 
effects and of the coupling of bound and continuum states 



A. S=0 pairing 

We shall carry out the calculation making use of a Gogny force [Q. In this case, the 
strongest part of the interaction v is the attractive S=0 component. The repulsive S=l 
component is much weaker and, within the present context, we shall ignore it. The spatial 
part of the S=0 component of the two-particle wavefunction ([^) is 

ipnn'ij{ri,r2) =< fi,f2\nn'lj] J = >s=o= 4>nij{ri)4>n'ij{r2) Pi{cos{e2i)), (6) 

where 621 is the angle between fi and r2. Inserting this expression into the (ri,r2) represen- 
tation of the function defined in Eq. we obtain the pairing density 

$5=o(ri, rl) = E (E U'niX'h) Mri)<Pn'iAr2) Pi{cos{02i)). (7) 

nn'lj q °^ 

This quantity is symmetric in the positions of the two particles, whereas the spin part, 
15* = >, is antisymmetric. Equation (|^) describes a correlated two-particle state, the 
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so-called Cooper pair wavefunction. Because of the f/, V factors, this wavefunction receives 
contributions essentially only from single-particle states around the Fermi surface. 

The pairing field is defined as product A(ri, = —v{\fi — r2|)$s=o(n_, ^2)- The pairing 
matrix elements are calculated, (cf. Eq. (^)), by projecting this field onto the two-particle 
wavefunctions (|^). For later discussion we introduce the Fourier transforms of the pairing 
field A and of the two-particle wavefunctions defined in Eq.(|]), 

A{k, R) = j rfVi2 A(fi, rs) e-*'^"^i^ (8) 

iljnn'ijik.R) = j ^/'nn'^i(rl,r2)e-^'=■'*^^ (9) 

where ri2 is the relative distance of the two neutrons and R is the center of mass of the pair. 
The pairing matrix elements can then be written 

Ann'lj = I d'k J d'R A{k, R) iJnn'uik, R) . (lO) 



B. Local Density Approximation 

Setting the nuclear potential U equal to zero, and allowing the radius Rb of the Wigner- 
Seitz cell go to infinity for a given value of the Fermi energy ep, one recovers the re- 
sults valid for uniform neutron matter. In fact, for R^ —>■ 00 one can substitute (pnij by 
{2k'^ / RbY^'^ji{kr), and the sums over n,n' by {Rb/nY J dk J dk' . In this limit, the matrix 
element of the gap defined in Eq.(^) becomes diagonal in k and A;', and the HFB equa- 
tions reduce to those of the BCS theory. The pairing density of Eq. (5) has only one spin 
component in this case, namely the S" = 0, and it is given by 

$™./(ri2) = / ^^^^fcVfce^^'-^-. (11) 

As expected, the Cooper-pair wavefunction depends, in this limit, only on the relative dis- 
tance ri2 of the two neutrons. This is also the case for the associated pairing field 
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Aumfiru) = -v{n2) J (12) 
The Fourier transform of this quantity is 

A.M^-l^,^^^^.irk-kv. (13) 

In this expression we have exphcitly indicated the dependence of the pairing gap on the 
Fermi energy, and the BCS relation UkVk = A{k)/2Ek has been used, the quasiparticle 
energy being 



E, = ^{h'ky2m - epY + A^^^/A:, ep). (14) 
The LDA to the pairing gap is obtained (cf. refs. P,p!0[|; cf. also |lT|) by replacing ep 



in Eqs. (p!3| ) and (|T^) by the local Fermi energy ep{R) = ep — U{R). The corresponding 
local Fermi momentum is given by h'^kl{R) = 2mep{R). In this way one obtains a space- 
dependent pairing gap AipiAik, R) = Aunif{k,ep{R)) and a space dependent quasiparticle 
energy, 



E,{R) = ^{h'ky2m - ep{R)y + A^/A;, ep{R)). (15) 



First- and second-order corrections in powers of h to the LDA have also been studied 12 . 



III. COMPARISON OF HFB AND LDA RESULTS 

In the following we present calculations of the quantum and LDA gaps A{k, R) for a 
Wigner-Seitz cell of 29 fm, and a Fermi energy ep = 13.5 MeV corresponding to 5o°°Sn and 
typical of the situation encountered at the density of about 3 xlO^^ g cm ( or p ^ 0.02 n 
fm~^) in the "inner crust" of a neutron star. As for the interaction v, we adopt the Gogny 
force in the ^5*0 channel p[. In the HFB calculations, we have performed the self-consistent 
calculation only for the pairing field, while instead of the Hartree-Fock field we have used a 
Woods-Saxon potential, without introducing an effective mass. The adopted Wood-Saxon 
potential (cf. Fig. 1) has a depth Vo = -31 MeV, a radius R = 7.5 fm and a diffusivity 
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a = 0.9 fm. These values are chosen in order to reproduce the neutron density in the cell, 
calculated in ref. [0. The resulting neutron density is equal to about Pint= 0.1 n fm~^ in 
the interior of the nucleus (R < 5 fm), and then drops to about pext = 0.018 n fm~'^ for 11 
fm < i? < 27 fm, until it goes to zero close to the boundary of the cell. We have included 
the single-particle levels up to 100 MeV, controlling the convergence of the results. We have 
also checked that in the absence of the Woods-Saxon potential, we obtain the pairing gap 
and specific heat of neutron matter at the same Fermi energy. The pairing gap in neutron 
matter Aunif{k,eF), calculated aX k = kp, is shown in Fig. 2 as a function of ep- In Fig. 2 
we also show the value of the local Fermi energy at three points in the Wigner cell. 

In neutron matter one can write ~ sin(/c^ri2)//cF'"i2 as one can see from Eq. (pUD 

taking into account the fact that the product f/^V^ is peaked at the Fermi energy, so that it 
oscillates with a wavelength given by A ~ 2Ti/kp. These oscillations are damped over a scale 
which is approximately determined by the coherence length ^ = h^kF/mirA ||T^ , p!5[| . The 



pairing gap calculated with the Gogny interaction in neutron matter at the Fermi energy ep 
=13.5 MeV is about 3.5 MeV (cf. Fig. 2), so that one obtain ^ = 3.2 fm. 

In Fig. 3 we show instead the square of the Cooper pair wavefunctions calculated with 
HFB in the presence of the nucleus at three values of R {R= 3 fm, 8 fm and 15 fm), 
as a function of r^. At large values of R, well outside the nucleus (cf. Fig. 3(c)), the 
wavefunction is essentially the same as calculated in uniform neutron matter at the same 
Fermi energy, ep = 13.5 MeV. The radius of the Cooper pair, given by the mean square 
relative distance, < >^^^, is about 4 fm, to be compared with the coherence length of 
3.2 MeV reported above. Inside the nucleus, at i? = 3 fm (cf. Fig. 3(a)), or on the surface, 
at R= 8 fm (cf. Fig. 3(b)), the wavefunction feels the presence of the nucleus, but its main 
features can still be understood in the same way, using the local values of the Fermi energy 
and of the pairing field. In this case the wavefunction shows some dependence also on the 
angle 6 between R and fu (especially for large values of ri2), and we show the results for 
6 = 45°. The local Fermi energy decreases going from inside to outside the nucleus (cf. Fig. 
2): it is equal to about 44.5 MeV at i? = 3 fm, instead of about 25 MeV at i? = 8 fm or 
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13.5 MeV at R= 15 fm. Therefore the local wavelength 2TT/kF{R) of the Cooper pair is 
much shorter inside than outside the nucleus, as can be seen comparing Fig. 3(a) with Fig. 
3(b) or Fig. 3(c). Moreover, the pairing gap Al/)a(^, ^^(-R)) calculated at the local Fermi 
momentum kpiR), denoted by A(/ci?(i?)), is much smaller inside the nucleus than on the 
surface or outside, (cf. Fig. 2 and Fig. 4(b)), so that the coherence length is larger at R = 
3 fm, because kp is larger, and A is smaller. Using the values of A{kp{R)) (1.6 MeV at R = 
3 fm and 3.5 MeV at i? = 8 fm, cf. Fig. 2), one obtains ^ = 12 fm at -R = 3 fm and ^ =4.1 
fm at -R = 8 fm. These values can be compared to the calculated r.m.s. values of 10.3 fm 
(at i? = 3 fm) and of 3.8 fm (at i? = 8 fm). 

The pairing gaps A{k,R) and A^DAik, R) are presented in Fig. 4(a) and Fig. 4(b) for 
different values of /c as a function of R. The value of A{k, R) depends on the angle between 
k and R, but we have found that this dependence is very small in the present case. The 
overall dependence on k and R is similar in the two calculations, but the HFB gap shows a 
smoother spatial variation than the LDA result. This is due to the rather large extension of 
the Cooper pairs, as discussed above and shown in Fig. 3, leading to important proximity 
effects, so that the difference in the gap between the interior of the nucleus and the outer 
part is smeared out. In the semiclassical case, we also show the values of the gap calculated 
at the local Fermi momentum as a function of R. It is seen that AiDAikpiR)) is equal to 
about 1.6 MeV and 3.5 MeV inside and outside the nucleus respectively, as already discussed 
above. 



IV. THE SPECIFIC HEAT 

A quantity playing an important role in the study of neutron stars, is the specific heat of 
the superfluid phase of the inner crust (cf. e.g. refs. 0, |T^, [|171)- This quantity is defined 

as 

_ l d<E> 

C.-y , (16) 
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where V is the volume of the Wigner-Seitz celL In the framework of our mean-field theory, 
the energy of the system is simply obtained summing over the quasiparticle states q: 

<E>=J2n,Eg, (17) 

the quantity = (1 + e^^'/-^^)"^ being the occupation number of the quasiparticle state q. 
The LDA result is obtained |]l6l taking the derivative of the quantity 



-^-^nkiR)Ek{R), (18) 

which is the semiclassical counterpart of Eq. (|l3). In this equation, the occupation factor 
nk{R) is given by nk{R) = (1 + e^'^''^^^/^))"^ where Ek{R) is given by Eq. (|^). 

In Fig. 5 we compare the specific heat calculated in HFB and in LDA. The semiclassical 
approximation grossly overestimates the specific heat in the presence of the nuclear potential 
at low temperatures, while it tends to the quantum result at high temperatures. 

In order to better understand the main source of this discrepancy, it is useful to write 
explicitly the expressions for at low temperatures, when A ^ T. It is also convenient to 
use the BCS approximation, which in the present case produces a value for the specific heat 
that is close to the HFB result. Using the fact that in BCS n = n', one can simply write 
Anij and ipnij instead of Ann'ij and 'ipnn'ij, and one obtains [|l^ 

= E A.,,e-(^-^--)^/^^^-^e-^-^/^. (19) 

nlj 

In the LDA we obtain instead for A ^ T, 

where we have used the relation given in Eq. (p!5|). 

For small temperatures, the dominant contributions to the semiclassical specific heat 
come from the regions of the phase space (fc, R) close to the line defined by the local Fermi 
momentum kpiR). As shown in Fig. 4(b), there is a clear difference in the values of 
A{kp{R)) inside and outside the nucleus, and the integral in Eq. (|^) roughly separates in 
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two contributions. Although the associated volume is small, the contribution from the inner 
region is very large, because the local values of the gap are small and are heavily weighted 
by the exponential factors. If one uses in Eq. pO| ) the HFB value of A{k, R) instead of the 
semiclassical values, one obtains an improvement compared the LDA, but the specific heat 
calculated in this way is still far from the quantum result. Considering instead the quantum 
expression ([ToD for the specific heat, and going back to Eq. ([ToD , we observe that the value 
of the gap for a given state is obtained averaging the value of A{k, R) over ipnijik, R). The 
main contribution to the specific heat comes from states close to Fermi energy, which are 
very delocalized, and are rather close to plane waves. On the contrary, Eq. ( ^Of ) uses the 
local value of the gap. This suggests to substitute the weighting factor e^^^^-^^'''^'^^'^ in Eq. 
(0) by e-<^^^A(fc,fl)>/T^ ^i^g^g 

LDAik' , R')'4>e{k' , R') (21) 

is the average of the semiclassical gap over the density associated with a single-particle state 
0e of energy close to e = ?i^k'^/2m + U{R). In the present case, it is sensible to approximate 
the wavefunctions close to the (positive) Fermi energy by plane waves e^'^^'^^'^'^' / ^/V , so that 

i^,{k\R')^^5{k' -k^iR')), (22) 

where the local momentum kL{R') is defined by 

^kl{R') + U{R') = 6, (23) 

and 

< ALDAik, R)>=l^J d^R'ALnAihiR'), R')- (24) 

Using the occupation factor q-<^lda>It (|20|), one obtains a specific heat which is 

very close to the exact result at all temperatures, as illustrated in Fig. 5. 

We conclude that the local density approximation to the pairing phenomenon over- 
estimates the spatial dependence of the gap, because it does not take into account the 
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non-locality of the Cooper pairs, and the associated proximity effects. Particular care has 
to be taken in estimating the thermal occupation factors, for which the LDA can lead to 
large errors at low temperatures. However, at least in the present case, it is possible to 
obtain an accurate approximation to the quantum HFB result by a simple averaging of the 
semiclassical expression of the pairing gap. 
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Figure captions 
Fig. 1 

The Woods-Saxon potential U (R) in the Wigner-Seitz cell of radius Rb = 29 fm is shown. 
We also indicate the Fermi energy ep =13.5 MeV, and the outer neutron density pext- 

Fig. 2 

The pairing gap Au„j/(A;, ei?) calculated with the Gogny interaction in uniform neutron 
matter at the Fermi momentum A;^ is shown as a function of the Fermi energy. Also indicated 
by the arrows are the values of the local Fermi energy ep^R) in the Wigner-Seitz cell at R= 
3 fm, 8 fm and 15 fm. 

Fig. 3 

The dependence of the square of the Cooper pair wavefunction $(-R, ri2) on the relative 
distance ri2 is shown for three different values of the center of mass of the Cooper pair in 
the Wigner-Seitz cell, = 3 fm (a), 8 fm (b) and 15 fm (c). 

Fig. 4 

The pairing gap A{k,R) is calculated as a function of R in the Wigner cell, using the 
HFB (a) or the LDA (b) equations. The gap remains constant for R> 12 fm. The various 
curves refer to different values of k, namely k — 0.25,0.75,1.25,1.75 and 2.5 fm~^ going from 
the top to the bottom curve. In figure (b), we also show the value of the gap calculated at 
local Fermi momentum kpiR) (dashed curve). 

Fig. 5 
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Specific heat of a Wigner cell containing a nucleus (in units of MeV fm calculated 
according to the HFB approximation, to the LDA approximation, and to the LDA corrected 
using the averaged pairing gap < ^lda > according to Eq. p^). 
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